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Abstract. We report the results of a search for the double-peaked blue-skewed infall signature in the bright-rimmed cloud core 
SFO 1 INE SMMl. Observations of the optically thick HCO+ and optically thin H''CO+ J=3-2 lines reveal that there is indeed 
a characteristic double-peaked line profile, but skewed to the red rather than the blue. Modelling of the dust continuum emission 
and line profiles show that the motions within SFO IINE SMMl are consistent with a collapsing central core surrounded by 
an expanding outer envelope. We show that the collapse is occurring at a similar rate to that expected onto a single solar-mass 
protostar and is unlikely to represent the large-scale collapse of gas onto the infrared cluster seen at the heart of SFO 1 INE 
SMMl. The outer envelope is expanding at a much greater rate than that expected for a photoevaporated flow from the cloud 
surface. The modelled expansion is consistent with the bulk cloud re-expansion phase predicted by radiative-driven implosion 
models of cometary clouds. 
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1. Introduction 

Spectroscopic evidence for collapse motions in low-mass star- 
forming regions has now become relatively widespread (e.g. 
Myers et al. 200(3 Evans I1999> . The classic spectroscopic 
signature of infall or collapse is an asymmetric blue-skewed 
double-peaked optically thick line profile with self-absorption 
at the systemic velocity (Myers et al. l2000> . The interpretation 
of this profile can be fraught with difficulty: unrelated clouds 
located along the line-of-sight can mimic the infall signature 
(Myers et al. I2000> . molecular outflows add extra complex- 
ity to the gas kinematics (Hogerheijde et al. '19981, chemical 
effects such as depletion may mask the infall signature com- 
pletely (RawUngs & Yates 120011 and rotation of the molecular 
core can completely reverse the blue-skewed profile to the red 
(Zhou [T995j . In order to conclusively demonstrate the presence 
of infall it is usually necessary to map the emission of both 
optically thick and thin lines across the molecular cloud core 
and compare their profiles to collapse models (e.g. Gregersen 
et al. EMOl Lee, Myers & Tafalla EMH . 

Despite these difficulties, collapse motions have been iden- 
tified in several protostellar and preprotostellar cores. The ob- 
jects that have been mainly investigated so far are Class & 
I protostars (Gregersen et al. 119971 Mardones et al. 1997i and 
preprotostellar cores (Lee, Myers & Tafalla [T9"99 2001 1. Most 
of these objects are located in nearby relatively isolated low- 
mass star-forming regions such as the Perseus, Serpens, Taurus 
and Ophiuchus molecular clouds (Mardones et al. .199?) . Few 
star-forming regions associated with the clustered star-forming 
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mode (Lada, Strom & Myers I1993I I have been searched for 
the presence of infall. Two notable exceptions include recent 
searches for collapse motions of gas associated with young stel- 
lar clusters themselves (Williams & Myers 1999; Williams & 
Garland 20Q2j and for infall within bright-rimmed clouds as- 
sociated with large HII regions (de Vries, Narayanan & Snell 
|2002J . 

Both types of search report limited success with collapse 
motions identified in two young clusters (Cepheus A and NGC 
2264) and one bright-rimmed cloud (SFO 18). It is the lat- 
ter type of object that we will focus on in this paper Bright- 
rimmed clouds are potential regions of triggered or induced star 
formation (Sugitani, Fukui & Ogura rr99Tl Sugitani & Ogura 
n'994j . whereby the collapse of the molecular gas comprising 
the cloud and the ensuing star formation process is initiated 
by the action of an external trigger In bright-rimmed clouds 
the external trigger is most likely the photoionisation-induced 
shocks that are driven into the clouds by the ionisation of their 
outer layers by nearby OB stars (e.g. Elmegreen 1991 ). These 
shocks compress the molecular gas of the clouds, triggering the 
collapse of the cloud and forming a dense, possibly quasistatic 
core at the cloud centre (Bertoldi 119891 Lefloch & Lazareff 
I1994lll995t . 

However, although de Vries, Narayanan & SneU (I2002> 
report the detection of an infall signature towards the bright- 
rimmed cloud SFO 18, it is the exception rather than the rule 
in their sample of bright-rimmed clouds. The overall blue ex- 
cess of the bright-rimmed cloud sample is less than that of 
the Class & I protostellar cores observed by Mardones et 
al. J1997t and Gregersen et al. (>2000i) . This phenomenon is 
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puzzling given the strong velocity gradients predicted by the 
Lefloch & Lazareff (1994i radiative-driven implosion model. 
However the small sample size of de Vries, Narayanan & Snell 
<2002< means that their result may not be statistically signif- 
icant. It is also possible that a flattened temperature gradient 
across the bright-rimmed clouds, caused by their external heat- 
ing, could mask the classic asymmetric infall signature (de 
Vries, Narayanan & Snell 2002i(. A wider infall survey of a 
larger number of bright-rimmed clouds, coupled with more re- 
alistic radiative transfer modelling of their temperature and ve- 
locity gradients is required to address the issue of infall in these 
clouds. Such a study would also permit the detailed investiga- 
tion of the kinematics of these clouds, which is important in the 
context of refining and expanding the existing radiative driven 
implosion models to take into account star formation processes 
(e.g. Bertoldi .1989,; Lefloch & Lazareff^ ,1994..1995j . 

We are currently undertaking such a survey of a statisti- 
cally significant sample of bright-rimmed clouds. As a fore- 
runner to the survey we present the results of a search for infall 
toward the bright-rimmed cloud SFO UNE. SFO UNE is lo- 
cated at the edge of the HII region IC 1848, approximately 4' 
NE of the bright-rimmed cloud SFO 1 1 (Thompson et al. 2004^ 
Sugitani, Ogura & Pickles 2002). SCUBA sub-mm continuum 
and JCMT CO observations reveal a dense core of molecular 
gas and dust located at the head of the cometary cloud (SFO 
IINE SMMl), whilst 2MASS images show that the dense 
dust core harbours a cluster of candidate young stellar objects 
and protostars (Thompson et al. 2004). The pressure balance 
and morphology of the cloud cloud suggest that the cloud is 
likely to be in the early collapse phase described by Lefloch & 
Lazareff' ( tl994j and is thus a good candidate in which to search 
for infaU. 

2. Observational procedure and results 

We obtained five spectra of the HCO+ J=3-2 Une at 267.5576 
GHz using the facility A-band heterodyne receiver at the James 
Clerk Maxwell Telescope (JCMT') during June 2002. The 
spectra were arranged in a 5-point cross-shaped grid cen- 
tred upon the sub-mm continuum peak of SFO IINE SMMl 
(q'(2000) = 02^ 51™ 53^;6, 5(2000) = -h60°07'00"). The FWHM 
beamwidth of the JCMT is ~22" at this frequency and the 
cross-shaped grid was spaced by 15". A single spectrum of the 
H'^CO+ J=3-2 line at 260.2555 GHz at the sub-mm continuum 
peak position was also obtained. The HCO^ line is usually opti- 
cally thick towards typical star-forming regions (e.g. Gregersen 
et al. 2000 ), whereas the less common isotopomeric line is usu- 
ally optically thin. 

All spectra were obtained using frequency-switched mode 
with a local oscillator switching offset of 8.2 MHz. backend 
spectrometer was set to its minimum bandwidth of 125 MHz to 
obtain the maximum possible velocity resolution of ~ 0. 1 km s 
per 94 kHz wide channel. The integration times for the HCO^ 
spectra and H'^CO^ spectrum were 20 minutes and 40 minutes 

' The JCMT is operated by the Joint Astronomy Centre on behalf 
of PPARC for the United Kingdom, the Netherlands Organisation of 
Scientific Research, and the National Research Council of Canada. 
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Fig. 1. Grid-map of the HCO^ J=3-2 spectra overlaid on an 
850 A/m SCUBA image of SFO IINE SMMl (Thompson et 
al. 2004 ). The map centre is at a(2000) = 02'^ 5 1™ 53!6, (5(2000) 
= -h60°07'00". 

respectively, which resulted in an rm.s. noise per 0.1 kms ' 
wide channel of 90 and 68 mK. 

The data were reduced using the Starlink package SPECX 
(Prestage et al. 2000). To improve the signal-to-noise ratio the 
spectra were binned to a velocity resolution of 0.2 kms"' . The 
HCO^ spectra are shown in Fig. [fl overlaid on the SCUBA 
850/vm image from Thompson et al. i2004t . The central HCO^ 
spectrum clearly displays a double-peaked profile in both the 
binned and unbinned spectra, although skewed to the red rather 
than the typical blue-excess infall profile (Myers et al. l2000> . 
A composite spectrum of the HCO^ and H'^CO^ spectra taken 
at the central position may be found in Fig.|2] The H'^CO^ line 
has a moderate asymmetry, with a clearly visible red wing. 

3. Modelling the dust continuum emission and gas 
kinematics 

The overall appearance of the HCO^ emission is puzzling. At 
the central position of the core the HCO^ emission is double- 
peaked, which is suggestive of infall, but the emission from 
the red peak is greater than that of the blue, perhaps implying 
that the core is expanding rather than collapsing. The opticafly 
thin H'^CO^ emission is singly peaked with a central velocity 
at that of the self-absorption dip, ruling out the possibility that 
the shape of the HCO^ emission is caused by two unrelated 
line-of-sight molecular clumps. Away from the central position 
the HCO^ spectra are more symmetric, without clearly visi- 
ble absorption dips. The outer spectra are also approximately 
the same brightness, although the southern and western spec- 
tra are slightly weaker than the remaining two spectra, which 
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Fig. 2. Composite spectrum of the HCO^ (thin) and the 
H'^CO^ (thick) J=3-2 lines observed toward the sub-mm con- 
tinuum peak of SFO 11 NE SMM 1 . Model spectra from the best 
fitting line model described in Sect.|3lare shown by dotted thin 
and thick lines. 

is perhaps due to the fact that the molecular cloud is sharply 
truncated in these directions (Thompson et al. '2004). The east 
and west spectra appear to be slightly velocity-shifted from one 
another, which may suggest that the core rotates about its north- 
south axis, although in these cases the emission toward the cen- 
tre of the core is usually symmetric rather than skewed to the 
red (Zhou[T995}. 

As mentioned in the Introduction, the interpretation of in- 
fall line profiles can be fraught with difficulty and in order 
to understand the kinematics of the molecular gas compris- 
ing SFO IINE SMMl we carried out detailed radiative trans- 
fer modelling of the sub-mm continuum and HCO^/H'^CO^ 
emission. As SFO 11 NE SMMl is only marginally elongated 
along the NW-SE axis in the SCUBA image and the line emis- 
sion suggests that the core is reasonably spherical with no 
strong velocity gradient we decided to model the emission 
using the RATRAN one-dimensional radiative transfer code 
(Hogerheijde & van der Tak 2000) . 

RATRAN solves for the hne and dust continuum emis- 
sion simultaneously. Firstly, we used the azimuthally averaged 
SCUBA 850 jum profile and the spectral energy distribution (at 
60, 450 and 850 jum) published by Thompson et al. (|2004| to fix 
the size, temperature and density profile of SFO 1 INE SMMl. 
In the modelling described in this paper we assumed the dust 
opacity model of Ossenkopf & Henning J1994t with thick ice 
mantles and 10^ years of grain growth. Our physical model for 
SFO 1 INE SMMl is that of a spherical core of radius 8 x 10'^ 
cm (0.26 pc) comprised of 10 logarithmically spaced radially 
symmetric shells with an inner cutoff at 5 x 10"' cm. Given 
that we possess information on the spatial distribution of the 
dust emission at only one wavelength (850 yum) we did not feel 
it was appropriate to model the density power-law index as an 
additional parameter and so we assumed a density profile r^^^^. 



We systematically explored a wide range of parameter 
space to find the best model fit to the data, including mod- 
els with inwardly-increasing, outwardly-increasing and flat 
(i.e. constant) temperature gradients. In all models we assumed 
that the gas and dust temperatures were equal, as is expected 
for high-density star-forming cores. The model flux profile was 
convolved by a 14" FWHM beam to match the SCUBA obser- 
vations. The best fitting model and the observed azimuthally 
averaged 850 yum flux profile are shown in Fig.|3l Fig. |3l also 
shows that there is an upturn in the 850 jim flux at large radii 
(>30"), which is more than likely due to the influence of the 
neighbouring core SFO 1 INE SMM2 lying ~50" to the south- 
east (Thompson et al. .2004i) . 

We could not reproduce either the observed spectral en- 
ergy distribution (SED) or the azimuthally averaged flux pro- 
file with a single-temperature model nor with an outwardly- 
increasing temperature gradient. An outwardly-increasing tem- 
perature gradient results in a much flatter flux profile than is 
observed, as the outer layers start to contribute more towards 
the total sub-mm flux. If the dust at the outer edge of the core 
is heated by the UV radiation from the nearby O-star it must 
comprise a relatively thin and low-mass layer that does not sub- 
stantially contribute to the 850 fim emission. By experimenta- 
tion we found that the best fits to the SED and flux profile were 
obtained by using a two temperature model with the inner three 
shells at 30 K and the remaining outer shells at 18 K. 

The best-fitting density profile for the two temperature 
model was found to follow the form riu-, (r) - 4.5 x 
10^ ('"max/'")^^^- Both the temperature and density compare well 
with the values derived from the single-component greybody 
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Fig. 3. The azimuthally averaged 850 //m flux profile and the 
best model fit to the profile. The 850 jim flux profile and its un- 
certainty is represented by the shaded area. The best model fit is 
shown by a solid line and the JCMT beam profile is represented 
by the dotted line. 
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model of Thompson et al. (I 2004> . The observed flux profile and 
SED are reproduced well by our two-temperature r^^^^ density 
profile model. The observed 60 and 450 yum integrated fluxes 
are reproduced to within 20% and the 850 /im integrated flux 
to within 10%. We estimated the sensitivity of the best-fitting 
model parameters to the data by progressively modifying the 
model parameters away from the best-fitting values. We find 
that the temperature (for each component) must lie within 2 K 
and that the outer density must lie within 3x10^ cm"^ of the 
best-fitting value. 

Using the dust continuum emission to fix the size of the 
core and its temperature and density profiles allows us to con- 
strain the model of the line emission very closely. The only 
free parameters in our subsequent modelling of the line emis- 
sion are the abundance of HCO^, the mean velocity of each 
shell and the velocity dispersion of the molecular gas. It is also 
possible to treat the ^^C/^^C abundance ratio as another free pa- 
rameter but in the modelling we carried out for this paper we 
assumed a standard '^C/'^C abundance ratio of 60 (Frerking, 
Langer & Wilson 1982 1. 

We modefled the HCO+ and H'^CO+ emission from SFO 
IINE SMMl along three broad hypotheses: pure infall, pure 
expansion, or a mixture of infall and expansion. The former 
set of models were included as a cross-check of our method 
as models comprised of pure infall should be unable to repro- 
duce the observed red-skewed profile. The red-skewed HCO^ 
line profile suggests the presence of expanding gas which may 
originate from either bulk expansion of the gas comprising 
the core or molecular outflows driven by protostars embed- 
ded within the core. There is some evidence that molecular 
outflows are associated with SFO IINE SMMl, by the pres- 
ence of moderate velocity CO line wings and a possible pho- 
toionised jet (Thompson et al.'2004l, but no other supporting 
evidence (e.g. spatially separated outflow lobes) has yet been 
found. Without spatially-resolved outflow lobes it is not pos- 
sible to model the gas expansion in terms of a molecular out- 
flow and so in all the expansion models described in this pa- 
per the expansion is modelled as a bulk expansion of the core. 
Nevertheless it is worth noting that any gas expansion inferred 
by the modelling may result from molecular outflows rather 
than a bulk expansion of the core. We will return to this point 
in Sect.m 

As expected, it was found that models comprised of pure 
infalling gas could not match the observed red-excess emission 
along the central line-of-sight nor the apparent symmetry of 
the outer lines-of-sight. Models consisting of a purely expand- 
ing core could not reproduce the level of asymmetry between 
the red and blue peaks of the central line-of-sight spectrum. In 
this latter case the blue peak is always much less bright than 
the red peak. For both sets of models we investigated whether 
velocity gradients across the model core (e.g. increasing infall 
velocity of each shell with decreasing shell radius or vice versa) 
could match the observed line profiles. In all cases the observed 
HCO^ emission could not be reproduced by either pure infall 
or pure expansion. 

We decided to restrict the models including a mixture of in- 
falling and expanding gas to a simple two-layer interpretation 
to simplify the modelling process and to reduce the number of 



free parameters in the model. This two-layer interpretation is 
similar to the simplified model described by Myers et al. J2000> 
but with spherical geometry rather than two uniform collapsing 
layers. For simplicity we will refer to the central portions of our 
spherical model as the "core" and the surrounding outer layers 
as the "envelope". We considered cases for both an expand- 
ing core combined with a collapsing envelope (e.g. NGC 2264 
IRSl; Williams & Garland 2002); and a collapsing core com- 
bined with an expanding envelope. In addition we modelled the 
subsets of these two cases with either static cores or envelopes. 

The best fit to the observed HCO+ and H'^CO+ spectra 
(i.e. with the lowest residuals) was produced by a model com- 
prised of a collapsing core surrounded by an expanding en- 
velope. Models consisting of an expanding core surrounded 
by a collapsing envelope produce a central blue-skewed pro- 
file and are not consistent with the observed spectra. The lim- 
ited number of free parameters in the model (HCO^ abun- 
dance profile, mean velocity and velocity dispersion) and the 
sharp discontinuity between the two layers meant that it was 
difficult to closely reproduce the depth of the HCO^ absorp- 
tion dip between the red and blue peaks of the central spec- 
trum. Nevertheless, the collapsing core and expanding enve- 
lope model closely matches the relative brightness of both 
peaks, their displacement from one another in velocity, the 
moderate red asymmetry of the H'^CO^ line and the relative 
symmetry of the outer spectra. 

We could not simultaneously fit any of the outlying HCO^ 
spectra in the grid-map with a single model and this probably 
reflects the limitations of assuming spherical symmetry in the 
model. However, the collapsing core and expanding envelope 
model produces the closest fit to the outer spectra of all mod- 
elled hypotheses. The models comprised of pure infall or pure 
expansion show pronounced asymmetry in their outer HCO^ 
spectra. Only by combining infalling and expanding gas mo- 
tions in the model can we reduce the asymmetry in the outer 
HCO^ model spectra to approach that seen in the observed 
outer HCO^ spectra. Also, given the close match of the model 
to the HCO'*" and H^^CO^ spectra along the central line-of- 
sight we are encouraged that our model interpretation is a good 
representation of the gas motions within the core. 

The largest discrepancy is in the spectrum to the north of 
SFO IINE SMMl, which is much broader and brighter than 
the model suggests. The eastern and western HCO^ model 
spectra are a much better match to the observations, although 
the western spectrum has a small velocity shift toward the blue. 
The difference between the eastern and western spectra may be 
due to the rotation of SFO IINE SMMl about its north-south 
axis, although the difference in the centroids of the two lines is 
too small (~ 0.2 kms"^) to show conclusive evidence for ro- 
tation. The model offers the best fit to the southern spectrum, 
reproducing the relatively broad red peak, the small possible 
self-absorption dip and the smaller blue peak visible in the ob- 
served spectrum. 

The best-fitting model specti-a of HCO+ and H^^CO+ to- 
wards the central line-of-sight of SFO 1 INE SMMl are shown 
in Fig. |2j whilst the model HCO^ spectra over the entire ob- 
served 5 -point grid are shown in Fig.|4] The velocity field that 
best reproduces the observed HCO ^ line profiles is of a central 
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Fig. 4. The HCO+ J=3-2 grid-map of SFO IINE (solid lines) 
and the best-fitting line model (dotted lines) described in Sect.|3] 



core collapsing at a constant velocity of 0.1 kms" surrounded 
by an outer envelope expanding at a constant velocity of 0.3 
kms"'. The central core has a radius of 0.1 pc, consisting of 
the inner eight model shells, and the outer envelope comprises 
the outer two shells of the model. 

It was found that the central HCO^ line profile is partic- 
ularly sensitive to the difference in velocity between the col- 
lapsing core and expanding envelope. Changing their relative 
velocities by only 0. 1 km s"' is sufficient to drastically alter the 
relative brightness of the red and blue peaks of the central spec- 
trum. The central HCO^ line profile is also sensitive to the rel- 
ative diameters of core and envelope. Good fits were obtained 
for an expanding envelope consisting of the outer three shells of 
the model, changing the envelope to comprise of the outer two 
layers or the outer four layers results in red peaks of the HCO^ 
central spectrum that are inconsistent with the observations. 

The HCO^ line profiles could not be reproduced by assum- 
ing a constant HCO^ abundance across the model and we found 
that smoothly increasing the abundance of HCO^ from 10"'' at 
the model edge to 5 x 10"** at the model centre gave the best 
results. The velocity dispersion in the model primarily affects 
the width of the resulting line profiles. We assumed a constant 
value of 0.3 kms"' for the velocity dispersion, which closely 
matches the observed linewidths of the HCO^ and H'-'CO^ 
lines. The model spectrum at the central position is narrower 
than the observed spectrum and this could indicate that the 
velocity dispersion increases towards the centre of the cloud. 
However, in the interests of restraining the number of free pa- 
rameters on the model we chose not to explore this possibility. 



4. A collapsing molecular core with an expanding 
outer envelope? 

In the previous section we described our physical model for 
SFO 1 INE SMMl as a dense central core of radius 0.1 pc col- 
lapsing at ~ 0.1 kms"', surrounded by an outer envelope ex- 
panding at a velocity of 0.3 kms"'. Here, we will explore the 
consequences of and offer a plausible explanation for the col- 
lapse and expansion of the molecular cloud core. 

Integrating over the entire model results in a total mass of 
25 Mo for SFO 1 INE SMMl, which is wi thin a factor of 2 of 
the greybody estimate of Thompson et al. ("2004). As the den- 
sity distribution in the model goes as r"^'^ and the model shells 
are logarithmically spaced, most of the mass lies in the outer 
shells of the model. For our collapsing core and expanding en- 
velope this translates into the core possessing roughly a third 
of the total mass (8 M©), whilst the envelope consists of the 
remaining two-thirds (17 Mq). 

The overall picture of SFO IINE SMMl gained from the 
modelling is of a low-mass central collapsing core surrounded 
by a relatively high-mass expanding envelope. It is unlikely 
that the double-peaked line profile is caused by two unrelated 
gas components at different velocities (e.g. the core and north- 
south ridge identified by Thompson et al. '20041 as the optically 
thin H'^CO^ line shows a single-peak at the systemic velocity. 
A possible explanation for the observed line profile is that the 
collapse motions result from infall onto one or more prototstars 
at the centre or the core, whilst the expanding outer layers rep- 
resent a photoevaporated flow caused by the illumination of the 
outer layers of the cloud by the nearby O-star. 

We may use the masses of each component, their radii and 
respective collapse and expansion velocities to estimate the 
mass infall and expansion rates to first order, allowing us to 
investigate this hypothesis. The simple relation M - Mv/r, 
relates the mass infall or expansion rate M to the mass of 
the core or envelope (M), the infall or expansion velocity (v) 
and the outer radius of the core or inner radius of the enve- 
lope (r). The infall rate Mi„ estimated from this relation is 
~ 7 X 10~^ Mq yr"' and the corresponding expansion rate Mexp 
is ~ 4 X 10"^ Moyr"^ 

The infall rate inferred from the model is consistent with a 
purely gravitational flow or free-fall collapse, i.e. Mi„ ^ cr" IG, 
where cr is the velocity dispersion (predicted to be 0.3 kms"' 
from the model) and G is the gravitational constant (Williams 
& Garland '20021. The infall rate is also approximately equal 
to that that estimated for a roughly solar-mass protostar (Reid, 
Pudritz & Wadsley 2002 Zhou et al. lT995l l. Is it possible tiiat 
we are observing infall onto a single protostellar object at the 
centre of SFO 1 INE SMMl? 

2MASS near-infrared images reveal a candidate cluster 
of class I protostars and YSOs at the heart of SFO IINE 
SMMl (Thompson et al. 2004) and so it is possible that we 
are observing gas that is either collapsing onto the cluster as a 
whole or perhaps onto individual protostars within the cluster. 
Hydrodynamical models of cluster formation suggest that as 
the cluster forms and contracts the steepening gravitational po- 
tential funnels an accreting gas flow onto the cluster (Bonnell 
et al. .200lJ . We note, however, that the mass infall rate within 
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SFO IINE SMMl is much lower than has been observed for 
other protostellar clusters (e.g. NGC 2264 IRSl; Williams & 
Garland 2002 1, where the infall rate can be up to 2 orders of 
magnitude greater than that for single protostars. For the gas to 
be undergoing a large scale collapse onto the cluster as a whole 
there must be a retarding mechanism in operation to slow the 
infall rate to that consistent with a single protostar. 

One of the predictions of radiative-driven implosion (RDI) 
models of bright-rimmed clouds is the presence of a photoe- 
vaporated flow from the surface of the cloud illuminated by 
flie nearby OB star (Bertoldi 1989, Lefloch & Lazareff' 1994 1. 
The expanding outer layers of the cloud, may in part be due 
to the presence of a photoevaporated flow from the surface of 
the cloud, although the large mass (2/3 of the total molecular 
cloud core mass) and size of the expanding layer are some- 
what difficult to reconcile with a photoevaporated flow alone. 
The mass loss rate from SFO 1 INE SMMl caused by the pho- 
toevaporated flow from the cloud surface was estimated by 
Thompson et al. i2004ij to be approximately 3.5 MoMyr"', 
which is roughly an order of magnitude lower than that im- 
plied by the mass and expansion velocity of the outer envelope 
(~4OM0Myr-i). 

The relatively high mass expansion rate impUes that the en- 
velope expansion is unlikely to be solely driven by a photoe- 
vaporative flow. There are two alternative hypotheses that may 
explain the presence of expanding gas in within SFO 11 NE 
SMMl: either molecular outflow(s) from protostars embedded 
within SFO 1 INE SMMl or a bulk re-expansion of the bright- 
rimmed cloud. The former possibility of outflow was men- 
tioned in the previous section: the expanding envelope of the 
model may be due to one or more molecular outflows mimick- 
ing the signature of the bulk expansion of the core. 

However, we do not see any evidence for the high-velocity 
line wings that are the usual indicators of molecular outflow 
in either the HCO^ spectra contained in this paper nor in the 
CO spectra of Thompson et al. f2004'l. Although the model 
of the cloud does not specifially consider a bipolar flow, we 
would expect the predicted mass of the expanding gas to agree 
with that contained in a bipolar outflow to at least first order 
The mass of outflowing gas predicted by the model (17 M©) 
is much greater than that typically observed in low mass out- 
flows, which is usually a few M© (e.g. Fukui et al. ll993l Richer 
et al. 1200 1. It is possible that the mass inferred by the model 
could be comprised of multiple overlapping flows, particularly 
in regard of the candidate cluster of protostars and YSOs found 
at the heart of SFO IINE SMMl. Although outflows cannot 
be completely ruled out as an underlying cause of the gas ex- 
pansion, the lack of supporting evidence for the presence of 
outflows in this bright-rimmed cloud (for example spatially re- 
solved outflow lobes) and the high outflow masses required to 
reproduce the model predictions substantially undermine the 
likelihood of this hypothesis. 

Another possibility is that the expanding gas is caused by 
a bulk re-expansion of the cloud. RDI models predict that 
this effect follows the maximum compression of the cloud by 
photoionisation-induced shocks. As the maximum compres- 
sion is reached, the cloud overshoots its equilibrium state and 
subsequently expands (Lefloch 8c Lazareff .1994j . Thereafter 



the cloud follows an oscillating and decaying cycle of com- 
pression and re-expansion until the gas reaches a quasi-static 
hydrodynamic equilibrium phase as a cometary cloud. Perhaps 
the expanding outer layers of the model cloud correspond to 
this expansion phase. 

Lefloch & Lazareff define a "radial momentum" for their 
RDI model clouds which is the integral of the gas momentum 
over the cloud volume. The radial momentum of our expanding 
envelope (0.8 M© pc Myr"') is well within the angular momen- 
tum limits evaluated by Lefloch & Lazareff ( 1994 1, even taking 
into account the relatively high mass fraction of the envelope to 
central core. Thus it seems that the expanding outer envelope 
is consistent with the re-expansion phase of the RDI models 
of bright-rimmed clouds, which is suggested to occure approx- 
imately .3x10^ years after the initial exposure of the cloud 
to ionising radiation (Lefloch & Lazareff [l994t . This value is 
consistent with the upper bound on the cloud age derived from 
the projected distance to the ionisation bound of the HII region 
(Thompson et al. l2004b . 

On balance, given the lack of supporting evidence for 
molecular outflows in SFO 1 INE SMMl we conclude that the 
most likely cause of the gas expansion is a genuine bulk ex- 
pansion of the cloud following the maximum compression of 
the cloud by photoionisation-induced shocks. The likelihood is 
that we are seeing SFO IINE SMMl in the beginning of its 
re-expansion stage, as this stage is the most consistent with the 
dynamic age of the associated HII region and the duration that 
the cloud has been exposed to the UV radiation from the nearby 
OB star. However, in order to confirm this hypothesis further 
observations aimed at confirming or denying the presence of 
bipolar molecular outflow are a high priority. 

5. Summary and conclusions 

We observed the bright-rimmed cloud core SFO 1 INE SMMl 
in sub-millimetre wave lines of HCO^ and H'^CO^ to search 
for the characteristic optically thick blue-peaked infall signa- 
ture (e.g. Myers et al. l2000t . The HCO+ emission toward SFO 
IINE SMMl is indeed double-peaked, but with an excess to- 
ward the red end of the spectrum rather than the blue. The 
H'^CO^ spectra show a single peak at the systemic velocity 
and this suggests that the double-peaked optically thick pro- 
file originates from a single gas core rather than a line-of-sight 
coincidence of two different velocity components. 

The line data were modelled with the radially symmet- 
ric radiative transfer model RATRAN (Hogerheijde & van der 
Ta kl2000t . SCUBA and IRAS HIRES data from Thompson et 
al. (2004"! were used to construct a detailed physical model of 
the temperature and density structure of the cloud core against 
which the Une emission could be modelled. The resulting phys- 
ical model of the source is a cloud core of radius 0.26 pc, with a 
density at the outer radius of 4.5 x lO"* cm"-* increasing toward 
the centre of the core as r"^^^. The SED and radial flux profile 
of the core suggest that the temperature of the core consists of 
two components: an inner warm region at 30 K surrounded by 
a colder envelope of 1 8 K. 

The HCO^ and H'^CO^ line profiles were found to be con- 
sistent with a combination of a central core coUapsing at a con- 
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stant velocity of 0.1 kms" , surrounded by an expanding en- 
velope moving outwards at a constant velocity of 0.3 kms"'. 
We derive the mass infall and expansion rates to first order 
by a simple consideration of the mass, velocity and radius of 
the collapsing and expanding regions. The mass infall rate is 
consistent with that estimated for a single solar-mass protostar 
(e.g. Zhou . 1995^ and it is possible that we are observing col- 
lapse onto a single member of the infrared cluster seen at the 
heart of SFO 1 INE SMMl (Thompson et al. [2004t . 

The expansion rate is an order of magnitude larger than can 
be sustained by a photoevaporative flow from the outer layers 
of the cloud (Lefloch & Lazareff 1994). The expansion is con- 
sistent with either the bulk re-expansion of the cloud predicted 
by the RDI models of Lefloch & Lazareff" ( I1994> or multiple 
unresolved bipolar molecular outflows. Given both the lack of 
supportive evidence for outflows and the relatively high mass of 
expanding gas predicted by the radiative transfer modelling we 
conclude that the likeliest scenario for the expanding gas is a 
bulk re-expansion of the bright-rimmed cloud. If so, SFO 1 INE 
SMMl is at the beginning of its re-expansion phase following 
the maximum compression of the cloud by photoionisation- 
induced shocks. 

Although our model of the cloud core as undergoing cen- 
tral collapse and outer expansion is consistent with the data, we 
stress that this is just one possible interpretation. Observations 
of higher sensitivity and spatial resolution are required to in- 
vestigate the kinematics of SFO 1 INE SMMl at smaller spatial 
scales and provide sufficient data for fully two-dimensional ra- 
diative transfer modelUng of the Une profiles (e.g. Hogerheijde 
& van der Tak l2000l PhiUips & Little l2000l l . In particular, 
the possibilities that the cloud core is rotating about its north- 
south axis and that the observed expansion may be in part due 
to bipolar molecular outflows must be investigated. Our forth- 
coming survey of a statistically significant number of BRCs 
will also reveal whether infall motions are as common within 
these clouds as more isolated low-mass star-forming regions. 
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